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ABSTRACT 

The energy spectrum of the lepton(s) in e + e _ — > tt — > ■ ■ ■ /£ + £~ • ■ ■ at next 
linear colliders (NLC) is studied for arbitrary longitudinal beam polarizations as a 
test of possible new physics in top-quark couplings. The most general non-standard 
form factors are assumed for yit, Zti and Wtb vertices to analyze new-physics 
effects in a mo del- independent way. Expected precision in determining these form 
factors is estimated applying the optimal-observable procedure to the spectrum. 
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1. Introduction 

The discovery of the top quark has completed the fermion spectrum required by 
the electroweak standard model (SM). It is still an open question, however, if 
the top-quark interactions obey the SM scheme or there exists any new-physics 
contribution. The top quark decays immediately after being produced JTJ since its 
huge mass m e t xp = 175.6 ± 5.5 GeV P] leads to a decay width r t much larger than 
/Iqcd- Therefore the decay process is not influenced by any fragmentation effects 
and the decay products carry lots of information on the top-quark properties. 

The energy distribution of the final lepton(s) in e + e" — > ti — > £ ■ ■ ■ /£ + £~ ■ ■ ■ 
turns out to be a useful tool to analyze the top-quark couplings |§. Indeed it 
has been frequently studied in the literature over the past several years in 
order to find observables sensitive to CP violation. To illustrate this point, it will 
be instructive to see how the spectrum is affected by non-conservation of CP in the 
production process: 

Since ti are produced mainly through j/Z exchange, their helicities would be only 
(-1 — ) or ( — h) if m t were much smaller than y^. Fortunately, however, this is not 

the case and we can expect copious (++) and ( ) productions as well even at 

\fs = 500 GeV.0 These states transform into each other under CP operation as 

CP | T =f) = | ± ±), which means that the difference N( ) — iV(++) could be 

a useful measure of CP violation 0-0. This important information cannot be 
drawn directly since the top decays too rapidly as mentioned, but is transferred to 
the final-lepton-energy distributions as follows: 

(1) The heavy top requires a large fraction (~ 70 %) of W bosons are longitudi- 
nally polarized in t — ► bW since Iry^^t ■ ~ mtbj^t when = ~ (e 
and k are the polarization and the four- momentum of W, respectively). 

(2) The produced b (b) is left-handed (right-handed) in the SM since m&/ 'y/s <Cl. 



S1 A rough estimate within the SM gives N(-+) : N(+-) : N( ) : N(++) is 5 : 3.5 : 1 : 1, 

where N(- ■ •) denotes the number of ti pairs with the indicated helicities 



(3) Because of (1) and (2), W +, s three-momentum prefers to be parallel (anti- 
parallel) to that of t(+)(t(— )), where £(• • •) expresses a top with the indicated 
helicity. Consequently £ + in the £(+) decay becomes more energetic than in 
the £(— ) decay, while it is just opposite for the t decay, i.e., £(— ) produces 
more energetic t~ than does. 

(4) Therefore, we expect larger number of energetic £ + (£~) for N( ) < 7V(++) 

(for N{ ) > N{++)). 

In realistic analyses, one should take into account that other source of non- 
SM effects may also exist. However, most of the above-mentioned articles focused 
on CP- violating effects in j/Zti vertices (production) only and did not assume 
the most general form for the interactions of ^ti, Zti and Wth. Therefore, in our 
previous paper || , we have performed a comprehensive analysis taking into account 
CP-violating and CP-conserving non-standard top-quark couplings contributing 
both to the production and decay process. 

In this paper, extending that work for arbitrary longitudinal e polarizations, 
we present a systematic way to determine the non-SM parameters describing the 



general 'y/Zti and Wth couplings. In our another recent paper |T0[ we have dis- 
cussed how the same process receives non-SM contributions from effective four- 
Fermi interactions. Therefore, with the present work we will complete a full anal- 
ysis of anomalous effects in top-quark interactions for polarized e + e~ beams in 
a model-independent way, where beyond-the-SM physics is parameterized by the 
577(3) x SU{2) x U(l) symmetric effective Lagrangian |12]| . 

This paper is organized as follows. First in sec. 2 we describe the basic frame- 
work of our analysis, and give the normalized single- and double-lepton-energy 
distributions. Then, in sec. 3, we estimate to what precision all the non-standard 
parameters can be measured using the optimal-observable method ||I3| . Adopting 
two sets of non-SM-parameter values we show in detail how effective the use of 
polarized beams could be for achieving better precision. Finally, we summarize our 
results in sec. 4. In the appendix we collect several functions used in the main text 
for completeness, though they could also be found in our previous papers [0, ||, [Tofl ■ 



2. The lepton-energy distributions 

In this section we briefly present our formalism, and then derive thereby the single- 
and double-lepton-energy distributions. 

We will treat all the fermions except the top quark as massless and adopt the 
technique developed in [|TJj]. This is a useful method to calculate distributions of 



final particles appearing in a production process of on-shell particles and their sub- 
sequent decays. This technique is applicable when the narrow-width approximation 
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can be adopted for the decaying intermediate particles. In fact, this is very well 
satisfied for both t and W since r t ~ 175(m t /miy) 3 MeV <C m t and TV — 2 GeV 
< M w . 

Adopting this method, one can derive the following formulas for the inclusive 
distributions of the single-lepton £ + and double-lepton £ + £~ in the reaction e + e~ — > 



tt 



e + e" £ + + • 



4 /■ .„ do- , , d 3 r. 



(e + e- - £ + r + 



d? P J{2f t )d?p'J{2 P y) 

= T*f d ^77T^ m ) o, (t ^ ^u) f% (t-+bl-u), (2) 

7 dQ t d 6 p e /(2p [ l) <Pp'J(2pl ) 

where i~£ and r t are the leptonic and total widths of unpolarized top respectively, 
and da(n, m)/df2t is obtained from the angular distribution of tt with spins s + and 
S- in e + e~ — > tt, do(s+, S-)/df2 t , by the following replacement: 

4 _> n ^ = + If 
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£ - m» = - [sT - ^ i^TPl. (3) 



(Exchanging the roles of s + and s_ and reversing the sign of n M , we get the single 
distribution of t~ .) 

In order to obtain the lepton spectra according to the above formulas we shall 
first calculate the tt-production cross section and their decay rates. 
ti production 

Let us start with the ti production. We can represent the most general ti 
couplings to the photon and Z boson as 

r vti = |«(P«) [r{A v + 5A v -(B v + 5B v )^}+ fo'J^ 1 * (6C v -5D v >y 5 )]v(pt), (4) 

where g denotes the SU(2) gauge coupling constant, v — 7, Z, and 

A^ = ^-sm9 w , P 7 = 0, A z = 1 (l-^sm 2 9 w ), B z = 1 . 
3 2 cos W v 3 ' 2 cos W 

Among the above form factors, 5A J>Z , 8B j Z , 8C 7t z and 5D 1;Z are parameterizing 
CP-conserving and CP-violating non-standard interactions, respectively. Note that 
we dropped two other independent terms proportional to {pt+PtY since their effects 
vanish in the limit of zero electron mass. 

On the other hand, interactions of initial e + e~ have been assumed untouched 
by non-standard interactions since their structures are well described within the 
SM: 



P; e+e - = -e^(p e+ ) 7 ^(Pe-), (5) 



7e + e vertex 



• Ze + e vertex 

r" Ze+ e- = J-&7rv(p e+ )^(v e + l5 )u( Pe -), (6) 
4 COS t>w 

where v e = — 1 + 4 sin 2 Ow- 

The angular distribution of polarized ti pair in presence of the above non- 
standard interactions is obtained after a tedious but straightforward calculation. 
The result is however a bit too lengthy, so we give the explicit form in the appendix 



and here instead we describe its structure rather qualitatively: 
First, the invariant amplitude can be expressed as 



M = Y,Cu jU IfM ( 7 ) 



u 



where 



4 = ^(p e+ )r>(p e -) (i — V, A) 



and 



r v,A,s,P = 1 ^ 7 ^ 75? q ^ g ^ 75 ( q = Pt - Pi ). 



Therefore consists of a number of terms whose coefficients are C* 7 Cj'j'. In 

the explicit formula in the appendix, we express C* z Cj// (1,1' = V,A), C*jCi>r 
(i ^ %' and /, V = V,A), CIjCvp (i, i', I = V,A) and QC,- S (i, i', I = V,A) as D, 
E, F and C respectively, and moreover we attach subscripts V, A and VA to D 
and E according to [I = V = V], [I = V = A] and [I = V, V = A], while F and G 
are classified by i = 1 ~ 4 according to their V/A structure.^] 
It is worth to notice that: 

• In the SM-limit only D v ,a,va and E V) a,va remain and all Fj = Gi = 0, 

• Non-zero Fj's are generated by the CP-violating form factors SD^ Zi 

• Contributions to Gj's are created by the CP-conserving form factors 5C lt z- 



For the initial beam-polarization we follow the convention by Tsai [IB]: 



P e - = +[N(e~, +1) - N(e~, -l)]/[iV( e - +1) + N(e~, -1)], (8) 
P e+ = -[N{e+, +1) - N(e+, -l)]/[N(e + , +1) + iV(e + , -1)], (9) 

where N(e~~( + \ h) is the number of e~(e + ) with helicity h in each beam.0 When 
the initial e~ and e + get polarized, and mix with each other since the 
spin (helicity) projection operator for u(p e -) and v(p e +) in the massless limit is 



* 2 More explicit formulas will appear in a separate paper |iq| . 

" 3 Note that P e + is defined with the opposite overall sign in some other papers (see, e.g., |]j 



(1 ± 75)/2. Then we obtain the cross section for arbitrarily-polarized e + e beams 
by replacing D v , D A , D VA , E v , E A , E VA , F { and (i = 1 ~ 4) with Dy\ D%\ 
D { *1 E { y\ E<x\ E { yl P« and Gf\ where 
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t and i decays 

We will adopt the following parameterization of the Wtb vertex suitable for the 
t — > W + b and t — >■ W~b decays: 
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«(p t ), (10) 
«(PS), (11) 



where Pl/# = (1 =F7s)/2, Vj& is the (£&) element of the Kobayashi-Maskawa matrix 
and k is the momentum of W. Because W is on shell,0 the two additional form 
factors were not taken into account. It is worth to mention that the above form 
factors satisfy the following relations ]T8|: 



ft R = ±fl 



L,R 



(12) 



where the upper (lower) signs are those for CP-conserving (-violating) contribu- 
tions.0 

Wlu couplings are treated within the SM as j/Ztt couplings: 
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(13) 
(14) 



" 4 Remember that we use the narrow-width approximation also for the W propagator. 
" 5 Assuming CP-conserving Kobayashi-Maskawa matrix. 



Assuming that f[ — 1, f% and f% are small and keeping only their linear 
terms, we obtain for the differential spectrum the following result: 



1 d 2 r e 1 + /3 W e 

-(t — > 01 V J = 



r t dxdcu y ' (3 W 



l + 2Re(/ 2 V^' ' 



u l + 2r 



(15) 



where x is the rescaled lepton-energy introduced in |JJ 



_2E t (1-P\ 1/2 
X ~ m t \l + [3) ' 

with Eg being the energy of I in e + e~ cm. frame, ui is defined as 



w = (Pt ~ Pe) /m 



2 /™2 
t ) 



i?^ is the leptonic branching ratio of t (~ 0.22 for I = e, /x), and 

= (l-r) 2 (l + 2r), r = (M w /m t ) 2 . 

An analogous formula holds for £ — > 6£~z/ with _/vf replaced by 
Lepton-energy distributions 

Now let us give the lepton-energy spectra in terms of x. Since we are going to ap- 



ply the method of optimal observables [13] in order to isolate various non-standard 
contributions, it is convenient to express the spectrum as a sum of known inde- 
pendent functions multiplied by coefficients parameterizing non-standard physics 
to be determined. In the following, we use as input data My/ = 80.43 GeV, 
M z = 91.1863 GeV, m t = 175.6 GeV, sin 2 6 W = 0.2315 § and ^ =500 GeV. 
1. Single distribution 

Adopting the formulas eqs. (p^, p4| , p~5D , keeping only linear terms in non-standard 
form factors, and integrating over Q t and the necessary top-quark-decay phase 
space, one obtains the following normalized single-lepton-energy spectrum: 

where er e e->tt = a tot (e + e~ — > tt) and ± corresponds to i^. The first term comes 
from the SM and the coefficients are 



the second term originates from the anomalous j/Ztt couplings (see eq.(|])) con- 
tributing to the production process 



4 = ax 5D { * ] - a 2 [8D$ - Re(G { * ] ) ] + a 3 Re(<5E>g) T £ W , 

and the third term comes from the non-SM Wtb couplings (see eqs. ([lO|JTT|) ) which 
influence the top-quark decay distribution (see eq. fllBD ) 



4 = Re(f 2 R ), c 3 - = Re(/ 2 L ). 

Here, 6D§\ VA are the non-SM parts of Dy \ VA , £W is a CP-violating parameter in 
the production process which is defined in a similar way as £ used in |5|, [7| : 

^W = 2He(F 1 w )oiS, 

and ai are defined as 

a x = - r/*)(3 - p 2 )a$l a 2 = 2 ^^41 «3 = 4 ag, 

with aH = l/[(3- (3 2 )D ( y'* ] + 2p 2 D$'* ) } (the superscript "(0)" denotes the SM- 
part) and 

(= 0.2074 in case of no beam polarization). On the other hand, the functions fi(x) 
are 

fi(x) = f(x)+ V ^g(x), f 2 (x)=g(x), f 3 (x) = 5f(x) + 77M 5g(x), (17) 

where f(x) and g(x) are functions introduced in ||, while 5f(x) and 8g(x) are 
functions derived in our previous work |7[| , which satisfy the following normalization 
conditions: 

■,(,>*-l and /^-/^-/«,^-a 

/(x) and describe the process with the standard top decays while 5f(x) and 
5(7 (x) come from the non-standard contribution to the decay process. Here let us 



remind readers that the cf term in (|T6|), which is proportional to g(x), originates in 
the spin dependent part of the lepton spectrum and would vanish if, for instance, 
hadronization effects would dilute the top-quark polarization. As explained in 
the introduction the lepton-energy spectrum should depend on the polarization of 
the parent top quark, that is the reason why all the non-standard effects in the 
production process manifest themselves as modification of the coefficient in front of 
g(x) for the normalized spectrum. We recapitulate these functions in the appendix. 

It should be emphasized that the coefficients cf contain both contributions from 
CP-conserving and CP-violating interactions, therefore their determination does 
not provide a direct test of CP invariance. However, as was discussed in ref. |7j one 
can easy combine measurements of the spectra for £ + and £~ in order to construct 
purely CP- violating observables like da + /dx — da~/dx. It is also worth to notice 
here that even though measurement of q does not disentangle CP-conserving and 
CP- violating interactions it allows for discrimination between non-standard effects 
originating from the production and those from the decay. 

The functions fi(x) are shown in fig.|Tj for unpolarized beams. Since /i,3(x) have 
P e ± dependence through r)(*\ we also present them in figsfj and |3] respectively for 
P e - = +1 vs P e + = 0/ + 1 and for P e - = -1 vs P e+ = 0/ - 1 as examples.0 



" 6 P e + = and +1(— 1) give the same rf*^ and consequently the same /i 3 (i) when P e _ = 
+!(-!)• 




Figure 2: The function f x (x) for P e - = +1 vs P e + = 0/ + 1 (solid line), for P e - = -1 
vs P e + = 0/ — 1 (dashed line) and for no polarization (dotted line). 




Figure 3: The function fz{x) for P e - = +1 vs P e + = 0/ + 1 (solid line), for P e - = —1 
vs P e + = 0/ — 1 (dashed line) and for no polarization (dotted line). 



2. Double distribution 

Applying the same algorithm as for the single spectrum one finds for the nor- 
malized double-lepton-energy spectrum the following formula: 

= *>■ < i8 > 

where the first term comes from the SM 

Ci = 1, 

the second and third terms are CP- violating non-SM contributions of ^y/Ztt and 
Wtb couplings respectively, 

C2 = &\ c 3 = Ue(f 2 R -f 2 L ), 

the fourth and fifth terms are both CP-conserving non-SM / y/Zti contributions 

c 4 = a[ 5D { * ] + a' 2 5D { * ] + a'^G^), 

c 5 = ai SDfi - a 2 [5D^ - Re(G'S* ) ) ] + a 3 Re(8D { *}), 
while the last term is CP-conserving non-SM Wtb contribution 

The corresponding functions are 

fi(x,x) = f(x)f(x) +77 w [f(x)g(x) +g(x)f(x) ] + rj' { * ) g(x)g(x) : 
f 2 (x,x) = f(x)g(x) - g(x)f(x), 
IJ ■•■■■'■) = 5f(x)f(x)-f(x)5f(x) 

+ V { * } [ $f(x)g(x) - f(x)Sg(x) + Sg(x)f(x) - g(x)Sf(x) ] 

+ r]'W[5g(x)g(x)-g(x)5g(x)}, 
f 4 (x,x) = g(x)g(x), 
f 5 (x,x) = f(x)g(x)+g(x)f(x), 
f 6 (x,x) = 5f(x)f(x) + f(x)5f(x) 

+ V { * ] [ ^f(x)g(x) + f(x)8g(x) + Sg(x)f(x) + g(x)Sf(x) } 

+ V '^[5g(x)g(x)+g(x)5g(x)\, (19) 



with r/M = p- 2 a { *l[ (1 + [3 2 )D { y'* ] + 2f3 2 Df* ] }{= 1.2720 for P e = P- e = 0) and 
being denned as 

a^[/T 2 (l + /?V(3-/3V*M4L 

a' 2 = 2(1 - /jyM) ag, a' 3 = 2(1 + /3V W ) 41 

/i,4,5,6(x, and f2,3(x,x) are respectively symmetric and antisymmetric in (x, x), 
which are signals of CP conservation and CP violation. Since f'4 and /s are both 
from the CP-conserving parts of the production process, we may recombine them, 
but we chose the above combination so that only f 5 remains in computing the 
single distributions. 

Here, as for the single spectrum, since for a given Cj there is no mixing between 
the production and decay processes, we will be able to judge if the non-standard 
contributions originate from the production or from the decay of top quarks. Fur- 
thermore, in contrast with the single spectrum, the coefficients q receive contri- 
butions either from CP-conserving (i = 1,4,5,6) or CP-violating (i = 2,3) inter- 
actions. Therefore determination of the coefficients provides a direct test of CP 
invariance. 

The functions fi(x, x) are plotted in fig.|] for unpolarized case. Since fi^^x, x) 
depend on P e ± through i]^ and/or we also show them in fig||for P e - = +1 vs 
P e + = 0/ + 1 (on the left side) and for P e - = — 1 vs P e + = 0/ — 1 (on the right side) 
as examples. It can be observed from the figures that the shapes of the functions 
/ 13 (x) and fx t 3fi(x,x) vary substantially with the polarization of the initial beams. 
Therefore it is justified to consider determination of the coefficients q through 
energy-spectrum measurements for various polarizations since one can hope that 
carefully-adjusted beam-polarization may increase precision of the analysis.^ 



" 7 Getting higher statistics is also a reason for considering polarized beams. 




Figure 4: The functions fi(x,x) defined in eq.(|19D for P e - = 



Pe + = 0. 




Figure 5: The functions fi,3fi(x, x) for P e - = +1 vs P e + = 0/ + 1 (on the left side) 
and for P e - = — 1 vs P e + = 0/ — 1 (on the right side). 



3. The optimal observables 

We are now ready to perform a numerical analysis, but let us first summarize the 
main points of the optimal-observable technique Suppose we have a cross 

section 

where /j(</>) are known functions of the location in final-state phase space and q's 
are model-dependent coefficients. The goal would be to determine q's. It can be 
done by using appropriate weighting functions such that / Wi(<p)S(<p)d(p = Ci. 
Generally, different choices for Wi(<p) are possible, but there is a unique choice so 
that the resultant statistical error is minimized. Such functions are given by 



w 



(0) = Ex u -/ J (0)/r(0), (20) 



where is the inverse matrix of My which is defined as 

M^jmm^. (21) 

When we take these weighting functions, the statistical uncertainty of (^-deter- 
mination through da /dip measurement becomes 



Ac, = y/X u a T /N , (22) 
where a? = f(da/d<p)d<j) and N is the total number of events. It is clear from the 



definition of the matrix My, eq.(^), that M^ has no inverse if the functions fi(4>) 
are linearly dependent, and then we cannot perform any meaningful analysis. One 
can see it more intuitively as follows: when fi{4>) = fj(4>) the splitting between q 
and Cj would be totally arbitrary and only c% + Cj could be determined. 
Numerical analysis 

We apply the above procedure to the normalized lepton-energy distributions 
derived in sec. 2, eqs( pJ^ , |T8|) . From the theoretical point of view, perfectly-polarized 
initial beams (P e + = P e - = ±1) are the most attractive. However, those are 
difficult to realize in practice, especially for the positron beam. We shall therefore 
discuss the following two cases: 



(1) p e+ = vs P e - = 0, ±0.5, ±0.8 and ±1, 

(2) P e+ = P e -(= P e ) = 0, ±0.5, ±0.8 and ±1. 

Before carrying out detailed computations, we shall briefly discuss how the 
statistical errors Act depend on P e ±. For this aim we have to check polarization 
effects in the lepton spectra. These spectra depend on P e ± through the coefficients 



Q and the functions /j in eqs. (P7|,[l9D as well, but the strongest dependence comes 



from the normalization factor since it is proportional to <7 eS -^tt which is 



<T e ^u ~ (3 - 1 ) (f + Pe-Pe+)Df - (P e - + P e+ )4 0) 



+ W 2 [ (1 + Pe-Pe+)Df - (P e - + P*)Ef ], (23) 

where Dy\ and Ey\ are the SM parts of Dy.A and Ey t A in eq.( p5|) . Neglecting the 
vector-type part of 'jee coupling v e (v e = — 1 ±4 sin 2 Q w is tiny for sin 2 Q w = 0.2315), 
we have 

DP = C(A 2 +A 2 z d' 2 ), Df = CB 2 z d' 2 , E® =2CA J A z d', Ef = 0. 
Since Ey > for sin 2 9y/ = 0.2315, negative polarizations increase a eS ^tt- 

The 

matrix M y - is proportional to u e e^tt through the normalization factor, which means 
that negative polarizations would reduce statistical errors, eq.(^2|), since the matrix 
Xij oc l/<7 ee _ > t£. As it has been mentioned, My- depends, to a certain extent, on 
P e ± also through Cj and 77 W in the functions /j, therefore even for nearly the same 
number of detected events (the same a eS ^tt) statistical errors may differ. However, 
the general tendency is consistent with this naive expectation as will be observed 
later in tables presenting our results. 
1. Single-distribution analysis 

First, we shall consider the single distribution. Using eq.(^) for da^jdx we 
can obtain Acf 3 , the statistical errors for the determination of function 
of the expected number of detected single-lepton events Ng. For a given integrated 
luminosity L and lepton-tagging efficiency eg one has = Bia eS _, t tL e eS , where 
L^ s = tiL (in fb _1 units) is the effective luminosity. In the following we use 



q = 0.6 and L = 100 fb 1 as an example of realistic experimental constraint^ and 
estimate a e e->tt within the SM by using a(s)(~ 1/126). 
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Table 1: Expected statistical errors in c^ 3 measurements and the number of the 
single-lepton- inclusive events Ne for beam polarization (1) P e + = vs P e - = 
0, ±0.5, ±0.8 and ±1, (2) P e+ = P e -(= P e ) = 0, ±0.5, ±0.8 and ±1 at 
\/s = 500 GeV. Ne has been estimated within the SM for ee = 0.6 and L = 100 
fb" 1 . 

In table |I| we present Act^ and Ne for the above ee and L with the described 
configurations of beam polarization. From table [IJ, readers might conclude that the 
use of polarized beam(s) is quite effective for providing higher precision, especially 
negatively-polarized beams seem to be most suitable since we have smaller Acts 
as anticipated in the above discussion. Indeed, this is the case for cf measurement. 
For instance, when Re(/^), Re(/,f) = ±0.1, then N$ D = \ct\/Act, statistical 
significances for an observation of ct, becomes 2.0 for P e - = — 1 and 3.3 for P e = 
— 1, which means we can expect 2a and 3a confidence level respectively. However, 
for a given set of non-standard couplings, the coefficients ct vary with polarization. 
Therefore we should discuss their N$d inevitably instead of statistical errors only. 

* 8 Assuming L = 100 fb" 1 is in fact quite conservative since the integrated luminosity as high 
as 500 fb -1 is being recently discussed as a realistic possibility in the context of the TESLA 
collider design for y/s = 500 GeV . 



We will consider the following two sets of the couplings (of the order of 15% of the 
SM strength) in tables g and |Q 

(a) Re(&4 7iZ ) = Re(SB 1:Z ) = Re(SC y>z ) = Re{8D^ z ) = 0.1, 

(b) Re(M 7 ) = Re(<5B 7 ) = Re(5C 7 ) = Re(fiD 7 ) = 0.1, 
Re(8A z ) = Re(8B z ) = Re(8C z ) = Re(5D z ) = -0.1. 
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+0.8 
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-0.5 


-0.8 


-1.0 


el 


0.39 


0.36 


0.28 


0.17 


0.38 


0.36 


0.34 


C2 


0.14 


0.16 


0.12 


0.05 


0.09 


0.06 


0.03 




3.03 


2.31 


2.25 


1.83 


4.10 


4.65 


4.96 


\<Z\/Ac£ 


1.11 


1.04 


1.01 


0.58 


1.00 


0.75 


0.49 


(2)Pe 
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+0.8 


+1.0 


-0.5 


-0.8 


-1.0 
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0.39 


0.28 


0.19 


0.17 


0.36 


0.35 


0.34 




0.14 


0.12 


0.06 


0.05 


0.06 


0.04 


0.03 


\4\/Acf 


3.03 


2.52 


2.47 


2.59 


5.20 


6.30 


7.01 


\C2\/M 


1.11 


1.13 


0.86 


0.81 


0.83 


0.68 


0.70 



Table 2: Statistical significance of cf measurement for beam polarization (1) P e + = 
vs P e - = 0, ±0.5, ±0.8 and ±1, and (2) P e+ = P e -(= P e ) = 0, ±0.5, ±0.8 
and ±1, and the parameter set (a) Re(&4 7 ) = Re(5A z ) = Re(5B 7 ) = Re(8B z ) = 
Re(5C 7 ) = Re(5C z ) = Re(5D 7 ) = Re(<SD z ) = 0.1 at y/s = 500 GeV. 



" 9 One may notice that certain entries (some of Cj coefficients) in tables g and [| are identical. 
Indeed two polarization scenarios considered here provide for these cases exactly same values for 
Cj. Therefore, comparing statistical significances for them one can see the net effect of different 
statistics, as the expected number of events is different for the cases. The same will also apply 
to tables || and |[ 
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0.07 
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0.11 
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0.01 
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1.97 


3.70 
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1.15 


1.07 


1.02 
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0.03 


0.24 


0.86 
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0.06 


0.02 


0.10 


(2)Pe 





+0.5 
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-0.8 
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0.17 


0.46 


0.59 


0.61 


0.08 


0.07 


0.07 




-4-10" 3 


0.11 


0.18 


0.19 


10~ 3 


0.01 


0.01 


\ci\/M 


1.33 


4.14 


7.86 


9.38 


1.20 


1.31 


1.44 




0.03 


0.97 


2.40 


2.95 


0.02 


0.12 


0.14 



Table 3: Statistical significance of cf measurement for beam polarization (1) P e + = 
vs P e - = 0, ±0.5, ±0.8 and ±1, and (2) P e + = P e -(= P e ) = 0, ±0.5, ±0.8 and 
±1, and the parameter set (b) Re(<£4 7 ) = — Re(M^) = Re(5B 7 ) = — Re(dB^) = 
Re(5C 7 ) = -Re(5C z ) = Re(<5D 7 ) = -Re(5D z ) = 0.1 at y/s = 500 GeV. 

These tables show that the use of negatively-polarized beam(s) is not always 
optimal: for the parameter set (a) a good precision in measurement is realized 
when P e < 0, but even in this case the precision in measurement becomes 
better for P e > or even P e = (table Q). Moreover in case (b) both and 
C2 get the highest precision for P e = +1 (table |3|). Therefore one should carefully 
adjust optimal polarization to test any given model of physics beyond the SM. One 
can conclude (as far as the coefficient sets discussed here are concerned) that the 
appropriate beam polarization for the set (a) provides measurements of at 5.0<r 
and 7.0a level for P e - = —1.0 and P e = —1.0, respectively. For the set (b) maximal 
statistical significance for c~2 determination is 6.6 and 9.4 for P e - = +1.0 and 
P e = +1.0, respectively. Since c\ it is seen that the maximal statistical 

significance for is much lower: 1.1 for the set (a) and 3.0 for the set (b). 
2. Double-distribution analysis 

We can perform similar computations for the double-lepton distribution. Re- 
sults are presented in tables f|, H an d We find again in table £| that negative 



polarizations give smaller As a result, \c3fi\/Acz } $ can be easily estimated 

from this table once Re(/^) and Re(/^) are fixed. On the other hand, 02,4,5 have 
polarization dependence themselves, so we need tables f| and ^ in order to draw 
a meaningful conclusion, where the statistical significance for 02,4,5 nas been pre- 
sented. Again some of Cj in tables [5] and |6| are identical as in the case of the single 
lepton channel. 
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0.31 


0.34 
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0.26 


0.22 


0.20 
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0.22 


0.19 


0.14 


0.12 


0.13 


0.10 


0.09 


Ac 6 


0.14 


0.12 


0.09 


0.08 


0.08 


0.07 


0.06 


N u 


1013 


893 


1063 


1278 


1641 


2260 


2775 



Table 4: Expected statistical errors in 02,3,4,5,6 measurements and the expected 
observed numbers of the double-lepton-inclusive events Nu for beam polarization 
(1) P e + = vs P e - = 0, ±0.5, ±0.8 and ±1, (2) P e+ = P e -(= P e ) = 0, ±0.5, ±0.8 
and ±1 at y/s = 500 GeV. Nu has been estimated within the SM for ee = 0.6 and 
L = 100 fb -1 . 

Among the coefficients for the double-leptonic spectrum, 02,3 are CP-violating 
parameters. Since C3 does not depend on the beam polarization as already men- 
tioned, one can just say (from table f|) that 3a effects could be observed for 
P e = —1.0 if Re(/ 2 R - fk)l^ > 0-18. On c 2 one has to conclude from tables § 
and ^ that for both sets of non-standard couplings its determination would not be 



easy for the assumed luminosity, as its statistical significance reaches at most 1.7. 
This is due to the smaller number of detected events in this channel as it could 
have been anticipated. Still we can say that the use of polarized beams is very 
helpful to increase precision. Indeed, if we are able to achieve L = 500 fb _1 as 
discussed in fll9"l , then |c 2 |/Z\c 2 would reach 3.8 for P e = — 1 in case (a) (the same 
value could be obtained for P e = +1 in case (b)), while we have only |c 2 |/Z\c 2 =1.4 
if the beams were unpolarized. 
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Table 5: Statistical significance of C2,4,s measurement for beam polarization (1) 
p e+ = vs P e - = ±0.5, ±0.8 and ±1, and (2) P e+ = P e -(= P e ) = ±0.5, ±0.8 
and ±1, and the parameter set (a) Re(<5A 7 ) = Re(5Az) = Re(5B 7 ) = Re(5P>z) = 
Re(5C 7 ) = Re(5C z ) = Re(5D 7 ) = Re(5D z ) = 0.1 at yfs = 500 GeV. 

045^ are CP-conserving coefficients. Concerning cq, 3cx-level measurement is 
possible for P e = —1.0 when Re(/ 2 R + f^)/^ > 0.18. On c 4 we are also led to a 
similar conclusion to c 2 , but C5 determination is different. That is, the statistical 



significance for C5 measurement can reach 2.0 for P e = — 1 (case (a)) and 3.3 for 
P e = +1 (case (b)). This is quite in contrast with that for c 4 , which is less than 2 
as one can see from tables |5] and |6|. 
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Table 6: Statistical significance of 02,4,5 measurement for beam polarization (1) 
P e+ = vs P e - = ±0.5, ±0.8 and ±1, and (2) P e+ = P e -(= P e ) = ±0.5, ±0.8 and 
±1, and the parameter set (b) Re(M 7 ) = — Re(M^) = Re(5B 7 ) = — Re(dB^) = 
Re(5C 7 ) = -R e (5C z ) = Re(SD 7 ) = -Re(<5D z ) = 0.1 at y/s = 500 GeV. 

4. Summary and comments 

Next-generation linear colliders of e + e~, NLC, will provide a cleanest environment 
for studying top-quark interactions. There, we shall be able to perform detailed 
tests of the top-quark couplings to the vector bosons and either confirm the SM 
simple generation-repetition pattern or discover some non-standard interactions. In 
this paper, assuming the most general (CP- violating and CP-conserving) couplings 
for jtt, Zti and Wtb, we have studied in a model-independent way the single- and 



the double-leptonic spectra for arbitrary longitudinal beam polarizations. Then, 
the optimal-observable technique [TB| has been adopted to determine non-standard 
couplings through measurements of these spectra. 

The method applied here, the optimal observables, allows to disentangle various 
non-standard contributions to the production process (e + e~ — ► ti) and to the decay 
(t — > Wb). Using the single- leptonic-energy spectrum for I and assuming non- 
standard couplings of the order of 15% of the SM strength, we have found that 
an appropriate selection of the initial-beam polarization may provide observable 
effects for non-standard corrections to the production process, \c\\jAc\, even at 
9.4a level when both e~ and e + beams are polarized and at 6.6a when only e~ 
beam is polarized. On the other hand, from the same spectrum measurement 
one can expect on non-standard contributions to the top-quark decay the statistical 
significance of the signal Nsd — \ c f\/Acf of the order of 3.0 and 2.0 for both beams 
polarized and only electron beam polarized, respectively. 

The direct application of the optimal method for the single spectrum does not 
allow for discrimination between CP-violating and CP-conserving non-standard 
interactions since their effects mix in coefficients of the spectrum, cf . However, as 
it was discussed in ref. one can easily combine measurements of the spectrum 
for £ + and t~ in order to construct purely CP-violating observables. 

In contrast with the single spectrum, utilizing the method of optimal observ- 
ables directly for the double-leptonic-energy spectrum one can separately determine 
and disentangle the CP- violating coupling from the production of ti pairs (02) and 
the one from the top-quark decay (03). For the typical strength of the non-standard 
couplings discussed here, the highest statistical significance for CP violation in the 
production and/or in the decay was estimated to be 1.7 for both beams polarized, 
while we found that the maximal signal from CP-conserving interactions in the pro- 
duction process (Ic^l/Ac^) could reach 3.3 and 2.3 for both and only electron beam 
polarized, respectively. For CP-conserving interactions in the decay the expected 
effect is lower, namely 1.6 for the statistical significance for both considered cases 
of maximal polarization. 



It should be emphasized that we have used in this study very conservative 
integrated luminosity, namely L = 100 fb _1 . That is, the luminosity considered 
now as realistic is by factor 5 larger. Therefore one may expect that even though 
we have not considered any background here and our analysis does not take into 
account any detector details (to a large extent they are not available yet), the results 
presented here should serve as a fair estimation of real signals for beyond-the-SM 
physics. 

To summary, we found (z) the use of longitudinal beams could be very effective 
in order to increase precision of the determination of non-SM couplings, however 
(ii) optimal polarization depends on the model of new physics under considera- 
tion, therefore polarization of the initial beams should be carefully adjusted for 
each tested model. For such optimal polarization the maximal non-standard sig- 
nal should be observable in the single-leptonic spectrum on the effects generated by 
contributions (both CP -conserving and CP -violating) to the production mechanism 
of ti pairs. On the other hand, the most challenging measurement would be the 
determination of CP-conserving contributions to the decay process. Since we have 
already carried out a similar analysis of possible consequences emerging from ef- 



fective four-Fermi interactions ee — > ti and t(i) — > bt r v(bt~ii') in [10], this paper 
completes a full analysis of modifications for lepton-energy distributions by non- 
standard interactions of the top quark in a model-independent way for polarized 
e + e _ experiments. 

The results presented here are the most precise ones which could be obtained 
from the single or double energy distribution alone. It will of course be possi- 
ble to achieve a higher precision by combing our results with other statistically- 
independent data. Among them, lepton angular distributions are very promising. 
Indeed what one could measure via the energy spectra are the real parts of the 
non-standard form factors, while we would be able to determine their imaginary 
parts by using, e.g., an up-down asymmetry to the top direction as shown in ||]. 
However, non-SM effects in the decay process were not taken into account in that 
study. The lepton angular distributions relative to the initial beam direction will 



also give us valuable information. Some analysis focusing on the CP violation in 
the production vertices has been made in [J20| . However, comprehensive analysis 
including non-standard effects both in the production and in the decay process for 
all measurable distributions of the ti decay products seems to be needed |H . 



Finally, let us give a brief comment on the effects of radiative corrections. All the 
non-standard couplings considered here may be generated at the multi-loop level 
within the SM. In fact, CP-violating couplings <5D 7i z and Re(/^ — f%) requires at 
least two loops of the SM, so they are negligible. However, CP-conserving couplings 
5A lt z, SB lt z, SC-y^z and Re(/.f + f£) could be generated already at the one-loop 
level approximation of QCD. Therefore, in order to disentangle non-SM interactions 
and the one-loop QCD effects it is important to calculate and subtract the QCD 
contributions from the lepton-energy spectrum, this is however beyond the scope 
of this paper. 
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Appendix 

The angular distribution of polarized ti pair is given by the following formula: 

( '" J 'e + e- - f(s+)f(s-)) 

D v [ {4m 2 s + (lq) 2 }(l - s+s_) + s 2 (l + s+s_) 



3f3a 2 
~ 16s 3 



+2s{ls + ls_ - Ps + Ps_) + 2lq{ls + Ps_ - Zs_ Ps + ) } 
+ D A [ {lqf{l + s + s_) - (4m 2 t s - s 2 )(l - s + s_) 

-2(s - Am 2 t )(ls + ls_ - Ps + Ps_) - 2 lq(ls+ Ps_ - /s_ Ps + ) 
-4 Re(D VA ) m t [ s(Ps + - Ps_) + lq(ls + + /s_) ] 
+2 lm(D VA ) [lqe(s+,s-,q,l) + ls-e(s+,P,q,l) + ls+e(s-,P,q,l)) 



+4 E v m t s(ls + + Is-) + 4E A m t lq(Ps + - Ps_) 
+4 Re{E VA ) [ 2m 2 t (ls+ Ps_ - Zs_ Ps + ) -Iqs] 
+4 lm(E VA ) m t [e{s + , P, q, I) + e(s_, P, g, Z) ] 

- Re(Pi) — [Zg s(ls+ - Is-) - {(Zg) 2 + Am 2 t s}(Ps + + Ps_) ] 
+2 Im(Fi) [se(s+,s_,P,g) + /g e(s+, s_, P, /) ] 

+2 Re(P 2 ) s(Ps+ Zs_ + Ps_ Zs+) 

- Im(P 2 ) — [ e(s+, P, g, Z) - e(s_, P, g, I) ] 

m t 

-2 Re(P 3 ) /g(Ps+ Zs_ + Ps_ Zs+) 

+ Im(P 3 ) e(s+, P, g, I) - e(s_, P, g, /) ] 
m t 

- Re(P 4 ) — [ Zg (Ps+ + Ps_) - (s - 4m 2 )(Zs + - Zs_) ] 
-2 Im(P 4 ) [ Ps + e(s_, P, g, Z) + Ps_e(s + , P, g, Z) ] 

+2 Re(Gi) [ {4m 2 s + (Zg) 2 - s 2 }(l - s + s_) - 2sPs + Ps_ 

+lq{ls + Ps- - Is- Ps+) ] 

-l m ( Gl ) ^-[e(s + ,P,q,l) + e(s-,P,q,l)\ 
m t 

- Re(G 2 ) — [ (s - 4m 2 )(Zs+ + Zs_) - Zg (Ps + - Ps_) ] 

-2 Im(G 2 ) [ Ps+e(s_, P, g, Z) - Ps_e(s + , P, g, Z) ] 

-Re(G 3 ) — [Zg (Ps+ - Ps_) - (s - 4m 2 )(Zs+ + Zs_) ] 
m t 

-2 Im(G 3 ) Zge(s + ,s_,g,Z) 

+2 Re(G 4 ) [ (s - 4m 2 )(Ps+ Zs_ - Ps_ Zs + ) + 2Zg Ps+Ps- ] 

+ Im(G 4 ) — (s - 4m 2 ) [ e(s+, P, g, Z) + e(s_, P,q,l)]], 
m t 

where (3(= \Jl — Am 2 / s) is the velocity of t in e + e~ cm. frame, 

P=Pe+ Pe (= Pt+Pt), l=Pe~Pe: Q = Pt ~ Pi , 

the symbol e(a,b,c,d) means e^ upa a^b u c p d a for eoi 23 = +1, 

D v = C[A 2 -2A J A z v e d' + A 2 z (l + v 2 )d' 2 + 2(A^-A z v e d')Re(5A 1 ) 
-2{A lVe d' -A z (l + v 2 e )d' 2 }Re(8A z ) }, 



D A = C \B\{\ + v 2 e )d' 2 - 2B z v e d'Re(SB 1 ) + 2B Z (1 + v 2 e )d ,2 Re(SB z )}, 
D VA = C [-A 7 B z v e d' + A Z B Z {\ + v 2 e )d' 2 - B z v e d\5A 1 )* 

+ (A 7 - v e d'A z )5B 1 + B z (l + v 2 )d' 2 (5A z )* 

-{A lVe d' -A z (l + v 2 )d' 2 }5B z ], 
E v = 2C[A 1 A z d' - A 2 z v e d' 2 + A z d'Re(5A 7 ) + (A 7 d' - 2A z v e d' 2 )Re(5A z ) ], 
E A = 2C[ -B\v e d' 2 + B z d'Re(5B 7 ) - 2B z v e d' 2 Re(6B z ) }, 
Eva = C[A 7 B z d' -2A z B z v e d' 2 + B z d'(5A 7 )* + A z d'5B 7 

-2B z v e d' 2 (6A z y + (A/ - 2A z v e d' 2 )5B z ], 
F x = C[-(A y -A z v e d')8D 7 + {A 7 v e d'-A z (l + vl)d' 2 }5D z ], 
F 2 = C[ -A z d'5D 7 - {A 7 d' - 2A z v e d' 2 )5D z ], 
F 3 = C[B z v e d'8D 7 - B z (l + v 2 e )d' 2 5D z ], 
F 4 ee C[-B z d'8D 7 + 2B z v e d' 2 5D z ], 

G 1 ee C[(A 7 -A z w e rf / )5C 7 -{v4 7 w e d , -A z (l+t; e 2 )d /2 }5C , z ], 
G 2 ee C [A z d'5C 7 + (A 7 d' - 2A^ e d ,2 )5C z ], 
G 3 ee C[-5^X5C 7 + J B z (l + t; 2 )d ,2 5C z ], 

G4 ee C [ B z d'5C 7 - 2B z v e d' 2 5C z } (25) 



C == 1/(4 sin 2 W ), d' ee s/ [4 sin^ cos 6 w {s - M|)]. 

In the above formulas, only linear terms in non-standard couplings have been kept. 
The functions f(x), g(x), 5f(x) and 8g{x) in eqs. (|TTD and (|T^) are given as 



and 




(for the interval I\, I4 



CAl-r) 2 , 
^i(l-x) 2 , 



(for the interval I2) 



(for the interval I3, Iq) 




(for the interval 1$) 



2 l + (3 . x(l + (3) 



g(x) = C 2 —rx + x — x In 



l-(3 ■ r(l-0) 

(for the interval ii, J 4 ) 
= C 2 { (1 - r + lnr)x + ^ (1 - r) 2 }, 

(for the interval J 2 ) 
= C 2 { (1 - x + lnx)x + 2(1 ^ } (1 - xf }, 

(for the interval I 3 , J 6 ) 

= C 2 x — m 



where 



1-13 1-13 

(for the interval 7 5 ) 



r 3 1 + _ 3 (1 + /3) 2 

^1 = 7TT77 — n — ) < - y 2 



2W (3 1 W (3 

and W = (1 — r) 2 (l + 2r) with r = (M w /m t ) 2 as defined in the main text, 

TO = C ? { \r(r + 8) - 2*(r + 2)i±| + ^(}±£) 2 



2 

r(l-/3) ' 



(for the interval I 1: J 4 ) 
= Cs{ ^(r - l)(r + 5) - (1 + 2r) lnr }, 

(for the interval J 2 ) 

= C 3 { ^{x - 1)(5 + 4r - 3x) - (1 + 2r) Ins }, 

(for the interval J 3 , J 6 ) 
^ r / x , 1 + 48x . 68 9 1 

(for the interval J 5 ) 
= C 3 [ 1 - (3 + 2(3 - (3)r + ^r 2 - |(1 - 2/?) (i±|) V 



+ {l + 2r + 2(l + /3)(r + 2Mln^-±|(], 



x{l + P) 
r(l-P) 
(for the interval ii, J 4 ) 



-(r - l)(r + 5) - (1 + 2r) lnr 
2 



+ (1 + /3)x{^(r - l)(5r + 1) - 2(r + 2) lnr} 

(for the interval I 2 ) 



2 
3 



4r - /3(2r + 1) + 2x{l -f3 + r{2 + (3)} 



+ ^(1 + 2f3)x 2 - {2r + 1 + 2(1 + f3){r + 2)x} \nx 



(for the interval I3, J 6 ) 
■(l + 2r)(^-ml±|) + ^ 



-2(r + 2)x{ T ^--(l + /3)ln^} 



1 + 







1-/5- 
(for the interval J 5 ) 



where 



C 3 = 



6 1 + /3 ^ 



W /3 1 + 2r 
The intervals ij (i = 1 ~ 6) of x are given by 

J i: r(l-/3)/(l + /?)<x<(l-/3)/(l + /3), 
I 2 : (l-(3)/(l + (3)<x<r, 
I 3 : r < x < 1, 

(/ li2)3 are for r > (1 - /3)/(l + /?)) 
h: r(l-(3)/(l + (3) <x <r, 
h: r<x<(l-P)/(l + P), 
J 6 : (1-P)/(1 + P)<x<l. 

(J 4 ,5,6 axeforr<(l-/3)/(l + /?)) 
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